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We compute the spectrum of the classical and quantum mechanical coarse-grained propagators for 
a piecewise linear discontinuous map. We analyze the quantum - classical correspondence and the 
evolution of the spectrum with increasing resolution. Our results are compared to the ones obtained 
for a mixed system. 
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In this work we consider systems with a highly complex regular dynamics exhibiting, as in the chaotic case, an 
endless hierarchy of phase space structures. Examples of such dynamics are found in polygonal billiards on a surface 
with spherical curvature l] and also in a class of maps introduced in [3, [3| that can be characterized as discontinuous 
piecewise linear. Typical phase space portraits of such maps are shown in Fig. ([T]). An infinite number of chains of 
elliptic islands pack the phase space. In the center of each island sits a stable periodic orbit of arbitrarily long period 
surrounded by a family of nested invariant curves which correspond to open orbits having an infinite but periodic 
sequence with the same period. The motion inside each island is regular. Their size decreases as the period increases 
and phase space takes a fractal structure. The boundary of the islands constitute the unstable set which is the infinite 
iteration of the discontinuity. The motion in these systems is non chaotic but highly irregular. In the case of the 
billiards this peculiar phase space is a consequence of the focusing mechanism on the sphere, while in the second case 
it is due to the discontinuity of the map. 

The corresponding quantum mechanical systems have been investigated in both examples. For the billiards it was 
shown Q that the quantum energy spectra follow a non universal intermediate statistics that can be understood 
in reference to the classical phase space plot, but no quantitative classical-quantum mechanical correspondence was 
derived. In Q the semiclassical regime of a piecewise linear map was explored using Gutzwiller's method of periodic 
orbit quantization. However, as pointed out by the authors, great difficulties arose in their semiclassical approximation 
when describing the eigenstates not supported by stable periodic orbits. The relative failure of Gutzwiller's approach 
was attributed to the ellipticity of the dynamics of the map. 

In the present contribution we attempt a different approach. We consider the dicontinuous mapping introduced in 
Q from the point of view of the spectral properties of the Liouville dynamics and study both the classical and the 
quantum Liouvillian at limited phase space resolution. The motivation is to compare the results to those well known 
that apply to chaotic or mixed systems 0, Q . 

In pure hyperbolic systems the asymptotic decay of classical correlation functions is exponential and the decay 
rates can be rigorously obtained from the Ruelle-PoUicott ( RP ) resonances ^]. To study the long time behavior of 
general chaotic and mixed systems, which are beyond the validity of the RP theorem, non rigorous methods have been 
developed to compute the resonances of the Frobenius-Perron(FP)propagator. All these approaches are based on a 
coarse grained Liouville dynamics of the density function in the limit of zero coarse graining [8] . In [9'] the blurring of 
phase space structures is implemented by adding a diffusive noise in the Liouville equation, which results in a coarse 
graining of the FP propagator and the limit of vanishing noise is finally considered . An alternative method developed 
in ^ uses a truncation of the infinite unitary FP operator to a finite dimension in a basis of functions ordered 
by increasing resolution. The eigenvalues of this non unitary NXN operator are calculated and the the limit A^ — > oo 
is taken. Since noisy propagation restores the quantum-classical correspondence, which is otherwise lost in short time 
scales, the computation of RP resonances for chaotic systems is a powerful tool to explore the link between classical 
and quantum mechanical dynamics and to look for the emergence of chaotic signatures in the quantum mechanical 
system. 

In the present work we compute the spectrum of the classical coarse grained FP and quantum ( Husimi) propagator 
for the discontinuous map on the sphere presented in Q , using the truncation method of [1, [1] . We follow closely the 
procedure that Haake and collaborators developed to study a system with mixed dynamics: the kicked top. We show 
that for low resolution classical and quantum mechanical propagators coincide, although the convergence is slow. We 
study the behavior of the eigenvalues and eigenfunctions of the truncated propagator with increasing resolution and 
compare our results to the ones obtained for a map with a mixed phase space [5| . 



The area preserving map considered in 3] acts on the angular momentum vector J = {Jx,Jy,Jz) = 
j (sin0 cos 0, sin sin(/>, cos0) of fixed length j. The corresponding phase space is the sphere, and cos 9 and 
the canonical variables. It has the form: 

M = R,{to) \Rx{fi)\ (1) 

that can also be rewritten as: 



M = R,{ij) Ry{-n/2) \R,{^i)\ Ry{n/2) (2) 

where \Ri{pi)\ denotes a rotation around the Ji axis, of angle ^ sign J^. Under the action of this map the points on the 
sphere with J^, < rotate in ^ while those with J^, > in — /i around the Jx axis. The whole sphere rotates then in 
uj around the Jz axis. That is, at each step, every point rotates linearly except for the ones on the great circle Jx = 
where the map becomes singular. The unstable set is constituted by all the images and preimages of this great circle. 
This is evidenced in the phase space portrait. The bounds of one island of each chain are tangent to the unstability 
line (/) = 7r/2,37r/2. 

The evolution of the phase space density p of the system is given by: 

Pn+l (cos e,4>) =V Pn (cos 9,4)) (3) 

where V is the FP propagator that is related to the invertible map ^ by: 

V p{cos9,(j)) = p{M-^ {cos 9, (f>)) (4) 

Eq.Q results from the formal integration of the Liouville equation: 

dt p^ C p (5) 

Defining the Liouvillians corresponding to each rotation in eq.Q, the FP operator corresponding to this map reads 
■P = exp[£jj^(^)] exp[/:;j^(_^/2)] exp[£|;j^(^)|] exp[/:;j^(^/2)] . 

The matrix elements of V can be computed in the basis of the spherical harmonics Yim{9, (f>), which span the space of 
phase space functions on the sphere and are ordered according to increasing resolution by index I. 
It can be easily shown that: 

ie^v[CR^(^)])irn,i'm' = 5i,i'Sm,Tn' exp[-imtj] (6) 



{eii-p[^\Rz{p)\])lni,l\-m' = 

d4>{( d cos 9 Y*^{9,ct^)Yi.^, {9, (p - p)+ I d cos 9 Y*^{9,ct^) Yv^, {9, ^ + 



{cos p for / = I' 

for even I - V (g) 

i/2 sinmp J -^d cos 9 Pi ^micos 9) Pi' _m' (cos 9) for odd / — /' 

where (ij„„/(ijj) are the Wigner's d-matrix and Pi^m{cos9) the Legendre functions. The last matrix element, which 
can be computed analytically, couples subspaces with different values of I. Therefore, a truncation of T-" up to a finite 
Imax leads to a non unitary matrix. 



The quantum version of the map in the Hilbert space of the wave functions spanned by the (2j + 1) eigenvectors 
of Jz, \im)^ is given by the Floquet operator: 



F = exp[-iwJ2] exp[-i/x| Ja;!] 

= exp[—iuiJz] exp[— i7r/2 Jj^] exp[— J2I] exp[i7r/2 Jj^] (9) 

where: 

\JA^Y.l=-,\m\\jm) 

The density operator p is now represented by the corresponding Husimi function Qp{6,(j)) — {jd(l)\p\ jOcj)) ( that is, 
its diagonal matrix element in the basis of coherent states) and its time evolution is given by: 



dtQp^Q Qp (10) 

The Husimi propagator exp[f/] is a unitary matrix of dimension (2j + 1)^. As for the classical case their matrix 
elements exp [5] = exp[C/j^^((^)] exp[t/j^ (-7r/2)] 6xp[5|fl^(^)|] e-xplGn (7r/2)] will be calculated in the basis of the spherical 
harmonics. 

As shown in [il] the Husimi propagator for rotations is identical to its classical counterpart. We are then left with the 
calculation of exp[f/|/j^(^)|])im,i'm', whose spectral representation is: 

eiL-p[Q\R,{p.)\]) ^ XI <3bmi)(im2| exp[-i(^|mi| - Imal)] P|j>„,)(j„2| (11) 

mi ,m2 = —j 

where Q and P ( the Q- and P-functions corresponding to \ jmi) {jm2 \ ) are respectively the right-hand and left- 
hand eigenfunctions of exp[t/|^^(^-)|]) with eigenvalue exp[— i/L((|mi| — |'7i2|)]. To express the matrix elements of the 
propagator ([Tl]) in the basis of the spherical harmonics, we need to calculate the scalar products < Yim\Q\jmi ){jm2 \ > 
and < Yijn\P\jmi >(ini2 1 > which are sums of Clebsch Gordan coefficients. We finally obtain: 

21 I rnin(j.,j+m) 

eMG\RAt^)\])irn,i'm' = J2 C^T™!-™).'™ ^ji^i~m)Mn. exp[-z/i( |mi | - |mi - m|)] (12) 

mi—rnax(~j-\-rn^—j) 

As in the classical case the matrix elements of the propagator are diagonal in m , but non diagonal in I. 

We can easily show that when exp[CJ] is truncated to a dimension N = {Imax + 1)^ with l„iax << 2j it approaches 
its classical counterpart In this limit the sum in eq. lfT^ can be approximated by: 



21 + 1 ■ 

exp[e|fl.(;.)|])im,i'm' »<5™,^,— -^[exp[i^m] ^ Cj™\_„)_;„Cj[;;\_„) ;,„-Hexp[-i^m] ^ (^]'ll,-^),iJ^]'^rL,^^),vJ 

Ta\ — —j mi — 1 

since |m| < Imax << 2j. Using then the asymptotic expression for / << j of the Clebsch Gordan coefficients |1Q |: 




CTnl.-n.Um^\l7:^^Yim{9,0) (14) 



( where 9 is the angle between j and the z-axis) and approximating jjfT Y^mi=-j /"i dcos9 we recover eq .([SJ. 
Therefore the matrix elements of the quantum propagator given by eq. (|12p coincide with their classical counterpart 
of eq .([5]) for low resolution an in this limit it will be equivalent to study the spectrum of the FP or the Husimi 
propagator. However the convergence of the quantum mechanical propagator to the FP propagator is much slower 
than the one found for the kicked top in 6], where the quantum mechanical corrections were of the order of 
l/2j + 1. This can be appreciated when comparing the quantal spectra of Fig. O with Imax ~ 20 and different 
values of j', with the corresponding classical spectrum of Fig|4]a). The convergence gets worse as Imax increases, as 
shown in Fig. [3]which displays quantal spectra for Imax = 40 and different values of j ( to be compared with Fig. |l]b)). 



We now diagonalize the truncated FP operator whose elements are given by eqs. (j6l7l8p in the subspace 

spanned by = {Imax + 1)^ spherical harmonics with < / < Imax- This corresponds to a resolution of phase 
space structures of area 4:Tr/({lmax + 1)^- Two sets of parameters uj, ^ are considered in the map eq.(IT]) with the 
corresponding phase portraits shown in Fig[TJ In the first case ( Case I, = /i = 7r(V5— 1)) few big elliptic islands are 
present, corresponding to orbits of low periodicity. In the second ( Case II, = 1, /i = 0.05 ) the quasi tori surviving 
in correspondence to the integrable case { fi — 0) are visible. 

In Fig. [H we display the eigenvalue spectrum of 'P^^'> corresponding to Case I for different values of Imax- For 
coarse resolution Imax = 20 most of the eigenvalues Xi concentrate inside a ring corresponding to 0.9 < |Ai| < 1. As 
the resolution increases, the outer ring becomes narrower ( for Imax = 70, 0.95 < |Ai| < 1) and the fraction of the 
eigenvalues located in the inner disk decreases. In Case II, the situation is very much the same for Imax = 70, that is, 
there is a very high concentration of eigenvalues inside a narrow ring and a few values lying in the inner disk. However 
for low resolution Imax = 20 the distribution looks different: most of the eigenvalues lie close to the unit circle. Other 
spectra have been computed for different values of the map's parameters lo and corresponding to different types 
of phase space portraits with more or less dominant islands. The characteristics of the spectra differ from case to 
case depending on the number of domains that can be resolved at each given resolution. This non universality was 
also remarked in [1] where the calculation of quantal spectra of triangular billiards on the sphere showed that the 
level distributions do not follow a universal statistics but can rather be understood in reference to the phase space 
portrait corresponding to each particular case. However the general features of spectra obtained with high resolution 
are the same in all cases. There is always a densely populated outer ring, corresponding to unimodular and almost 
unimodular eigenvalues, and a small fraction of the eigenvalues lying in the inner disk. The eigenvalues move about 
when increasing Imax and no 'frozen' eigenvalues can be individualized. 

A better understanding of the spectrum is achieved by plotting the distributions of the N eigenvalues as a function of 
their modulus for different resolutions. As shown in FigslH](for Case I and II, and Imax = 50, 60, 74) we can distinguish 
three regions: a peak at 0.95 < |Ai| < 1 , an intermediate domain corresponding to approximately 0.85 < \Xi\ < 0.95 
which looks very much the same for all resolutions and a small number of values with |Ai| < 0.85. While the peak 
corresponding to 0.95 < |Ai| < 1 grows for increasing values of Imax, the number of states with |Ai| < 0.95 remains 
fairly constant, so that the fraction of non unimodular eigenvalues decreases with increasing resolution. Moreover, 
the numerical results seem to imply in the intermediate region a relationship rt(|A|) ^ |A|" where the exponent a is 
independent of the resolution. At this stage we cannot interpret this observation. The eigenfunctions corresponding 
to these three regions are of different nature. Some examples are represented in FigOfor Case I and Imax — 70 and in 
FiglUfor Case II and Imax — 70. Eigenfunctions with unimodular eigenvalue are shown in Fig. [7^) and FigEK): they 
are supported by the elliptic islands that can be resolved at the corresponding Imax- Eigenfunctions corresponding 
to quasi unimodular eigenvalues ( 0.9 < |Ai| < 0.95) have the structure displayed in FiglTjD) and Figl5)3) . Their 
amplitude is evenly distributed over a large number of islands. Finally for |Ai| < 0.85 the eigenfunctions are of the 
type shown in Figl7t),d) and FiglHl:),d): they locate along the unstabiliy line and along curves which are their first 
iterations by the map M. 

On the basis of these numerical results we can describe the V^^'^ spectrum, consisting of unimodular, almost 
unimodular and non unimodular eigenvalues, and its evolution with N as follows. Unimodular eigenvalues correspond 
to eigenfunctions which locate in the elliptic islands visible at the given resolution: when N grows more and more 
elliptic structures are resolved and their number increases. The group of almost unimodular eigenvalues correspond 
to eigenfunctions which are not fully localized yet and spread over several islands. The fact that their number 
remains approximately constant with iV, showing a compensation between the migration of existing eigenvalues to 
the unit circle and the appearance of new ones, related to still unlocalized states, can be explained by the existence 
of an infinite hierarchy of phase space structures which are gradually resolved. Finally, we saw that eigenvalues with 
moduli I Ail < 0.85 correspond to eigenstates sharply localized on the unstability line. These states are present for 
any value of N. However, their number seems to be approximately constant, so that they will represent a vanishing 
fraction of the total spectrum in the limit N oo. 

As expected, the observed evolution of the 7^(^) spectrum with N is different from the one obtained for the integrable 
system (corresponding tho /i = in our map) , where all eigenvalues migrate to the unit circle in the N —^ oo limit 
and unitarity is recovered. The situation is also different from the one observed for chaotic systems, where 7^(^) has 
an essential spectrum inside a circle of radius (r > 0) and a point spectrum inside a disk r < |A| < 1, with some 
eigenvalues (that can be related to the RP resonances) persisting in their positions for increasing values of iV 
In the system under study, non unimodular eigenvalues exist and can be interpreted, as in the chaotic case, as an 
effective dissipation due to the loss of probability to finer unresolved structures. However 'freezing' which in the case 
of chaotic dynamics is attributed to the self similarity in phase space, is not observed. Thus the 'resonances' will not 
be a characteristic of the classical system but rather will depend on the resolution. In addition, the fraction of non 



unimodular eigenvalues vanishes in the limit N —^ oo, since it is related to the measure of the unstable set. 

Summarizing, we computed the matrix elements of the FP operator 7^(^) and of the quantum Husimi propagator 
for discontinous piecewise linear systems and showed that they coincide in the limit of low resolution. We investigated 
the evolution of the spectra of the truncated propagator with N. We observed significant differences with 

respect to the results obtained for systems with chaotic or mixed phase spaces, where this truncation procedure is 
used to get the RP resonances related to the decay rate of correlations. One is the absence in our case of stable 
eigenvalues persisting in their positions for iV ^ oc. Another fundamental difference is the fact that the fraction of 
the non unimodular eigenvalues, present for all values of the resolution and corresponding to eigenfunctions localized 
on the unstability line and their images , decreases for increasing values of the resolution. This is can be taken as an 
indication of the zero measure of the unstable set. 
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FIG. 1: Phase space portrait of map eq.([T} for a) lo = fi — 7r(\/5 — 1)) ( Case l),h) uj — l/i — 0.05( Case II). 
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FIG. 2: Quantal spectrum (Case I) corresponding to Imax ~ 20 for j=100,200,300. 
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FIG. 3: Quantal spectrum (Case I) corresponding to Imax = 40 for j = 100, 200, 700. 
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FIG. 4: Classical spectrum (Case I) for Imax = 20, 40, 60, 70. 
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FIG. 5: Classical spectrum (Case II) for Imax = 20, 40, 60, 70. 
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FIG. 6: Number of eigenvalues (in a logarithmic scale) as a function of ln\\i\ for different values of Imax in a) Case I and b) 
Case II. 




FIG. 7: Eigenfunctions of the classical propagator 'P^''' ( Case I , Imax = 70) corresponding to eigenvalues with a)|A| = 1, 
b)|A| = 0.97, c)|A| = 0.85, d)|A| = 0.72. 
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FIG. 8: Eigenfunctions of the classical propagator 'p''^' ( Case II , Imax = 70) corresponding to eigenvalues with a)|A| = 1, 
b)|Aj = 0.95, c)|Aj = 0.88, d)jA| = 0.75. 



* Electronic address: |Spina@ tandar. cne a.gov. m] 



saraceno@tandar. cn ea.gov. ar I 



Electronic address: 
[1] Spina M.E and Saraceno M. 1999 J.Phys.A 32 7803 
[2] Scott A.J. , Holmes C.A. and Milburn G.J. 2001 Physica D 155 34 
[3] Scott A.J. and Milburn G.J. 2001 J.Phys.A 34 7541 
[4] Spina M.E. and Saraceno M. 2001 J.Phys.A 34 2549 

[5] Weber J., Haake F. and Seba P. 2000 Phys. Rev. Lett. 85 3620 ; Weber J., Haake F., Braun P., Manderfeld C. and Seba 

P. 2001 J. Phys A 34 7195 
[6] Manderfeld C., Weber J. and Haake F. 2001 J. Phys A 34 9893 
[7] Ruelle D. 1986 Phys. Rev. Lett. 56 405 

[8] Fishman S. and Rahav S. 2002 Lectures Notes in Physics: Dynamics of Dissipation, ( eds.Garbaczewski P. and Olkiewicz 

R., Springer- Verlag, Berlin) 
[9] Garci'a-Mata I., Saraceno M. 2004 Phys. Rev. E 69 056211 
[10] Varshalovich D.A., Moskalev A.N. and Khersonskii V.K. 1988 Quantum Theory of Angular Momentum (World Scientific) 



